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SUMMARY 

Induction heating of a low -density cylindrical plasma is investigated experimentally 
and analytically; in particular, the ohmic power dissipation and energy storage in the 
plasma are studied. The experimental values for both quantities are in good agreement 
with the theoretical results. 

The theoretical treatment consists of a linearized analysis of the coupling of a sole- 
noidal radio -frequency magnetic field and a steady cylindrical plasma. The plasma is 
assumed to have uniform electron density, electron temperature, and resistivity. An 
important feature of this analysis is that the inertia term is retained in the Ohm's law of 
the plasma. The ratio of this term to the resistive term is of order wr (the product of 
the applied frequency and the electron collision time) . The effect of the inertia term is 
to permit deeper penetration of the applied field into the plasma. The power dissipated 
and average energy stored in the plasma over a cycle are calculated by evaluating the 
complex Poynting vector at the plasma surface. The effect of the energy storage is to 
reduce the inductance of the radio-frequency coil. An interesting result of this analysis 
is that for low pressure plasmas this inductance change is a function of electron density 
only. This result indicates that a plasma may possibly be generated and diagnosed with 
the same radio-frequency coil. 

The experiments were conducted in helium at pressures of 3 to 20 microns with a 
17. 5-megahertz, 2-kilowatt radio -frequency power supply. The discharge tube was 
4. 8 centimeters in diameter. The power was coupled to the plasma by a 4 turn 7. 7- 
centimeter-long, 7.5-centimeter-diameter Faraday shielded radio-frequency coil. The 
power dissipated, the average energy stored in the plasma, and the radio -frequency mag- 
netic field were determined from detailed electrical circuit measurements. To compare 
theory and experiment, line ratio spectroscopy and microwave interferometry were used 
to determine the electron temperature and density, respectively. Power transfer effi- 
ciencies as high as 73 percent were obtained with 650 watts being coupled to the plasma. 

♦ £h 

Preliminary results of this investigation were presented at the 8 annual APS 

Plasma Physics meeting, Nov. 2 to 5, 1966, Boston, Mass. 



INTRODUCTION 


The excitation of gases by means of an electrodeless induction discharge has a long 
history dating from the early work of Hittorf (ref. 1) who first observed that an electrical 
discharge could be produced in a gas tube by induction from a coil surrounding the tube. 
The radio -frequency (r-f) electrodeless discharge has been useful through the years as a 
spectroscopic light source (ref. 2). In recent years, interest has been renewed in the 
r-f induction heating discharge as a plasma or ion source for plasma propulsion devices 
and in connection with energy addition to flowing gases (ref. 3). The use of r-f induction 
heating as a plasma source is advantageous in that acceptable power transfer efficiencies 
are attainable while the problem of electrode erosion is eliminated. An analysis of the 
coupling between a solenoidal r-f magnetic field and a steady cylindrical plasma may 
also aid in understanding the basic mechanisms in low-pressure plasma acceleration de- 
vices such as the travelling magnetic -wave accelerator (ref. 4). 

The present theoretical and experimental analysis of r-f induction heating was con- 
ducted to determine the basic plasma properties and power transfer efficiencies as func- 
tions of the external circuit parameters. The theory is developed and solutions are ob- 
tained for the power coupled to the plasma and the inductance change of the r-f coil. The 
experimental arrangement and diagnostics needed to compare experiment with theory are 
then described. Finally, the experimental and theoretical results are compared and dis- 
cussed. 


THEORY 

Assumptions and General Equations 

In this treatment, the plasma and r-f coil are considered to be coaxial circular cyl- 
inders of sufficient length that end effects may be neglected. The coupling is considered 
to be purely inductive and, therefore, there are no axial electric fields or currents. 
Furthermore, it is assumed that only alternating currents appear with no charge buildup. 
For tube radii in the centimeter range and megacycle field frequencies, displacement 
currents are negligible and Maxwell's equations take the form 
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where B = 


m 0 h 


and D = e^E. 


(Rationalized SI units are used throughout this report, 


and all symbols are defined in appendix A. ) The dielectric constant Cq and magnetic 
permeability /Iq are assumed to be the same in the plasma as in free space. 

The linearized Ohm’s law for the plasma, including electron inertia, is 
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where the electron density, n and resistivity 77 = are assumed to be uniform 

e 2 

n e q 

throughout the plasma. In actuality, the electron density and resistivity are not uniform 
along the plasma radius. However, the assumption made is that these spatially dependent 
quantities may be replaced by a uniform electron density and resistivity which represents 
the effective values of these quantities (Eckert, ref. 5). 

In the subsequent analysis, only the azimuthal component of equation (5) is used. 

— ► >•**> 

The previous assumptions of no charge buildup and j = 0 require that j = 0. Thus, no 

z r 

loss in generality results from ignoring either the j x B or terms in the generalized 
Ohm's law (ref. 6), since these terms are identically zero in the azimuthal direction. 

The only other term in the generalized Ohm's law (ref. 6) that may not be of higher order 
and is neglected in equation (5) is the v x b term. This term will be present only in the 
presence of a direct -cur rent magnetic field. However, for sufficiently small direct- 
current magnetic fields this term may be neglected. 


Formation of Complex Poynting Vector 

The complex Poynting vector will now be formed in order to obtain solutions for the 
power dissipated in the plasma and the energy stored in the plasma. Since equations (1) 
to (5) are linear, the following separation of variables is made 


E(r, t) = E(r)e ia>t 

(6a) 

B(r, t) = B(r)e ia,t 

(6b) 

3 (r, t) = J(r)e lwt 

(6c) 


Using equations (6a) to 6(c) reduces equations (3), (4), and (5) to 
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where E(r), B(r), and j (r) are complex time independent vectors whose amplitudes cor- 
respond to the physical peak values. In most of the subsequent equations only their de- 
pendence on r will be understood. 

The complex conjugate of equation (8) is 

Vxb*=jU 0 T* (10) 

Forming the scalar product of equation (10) with E and equation (7) with -B and add- 
ing yields 




+ iwB • B 


* 


which becomes, on regrouping terms and using equation (9), 



Equation (11) is an extension of the result given by Abraham and Becker (ref. 7, p. 202, 
eq. (9)). * After integrating equation (11) over the volume and using the divergence 
theorem, the result can be written 


- J- fit. ExbI. 
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dS = W + 2iwU 


(12a) 


where 


A The authors are indebted to Dr. F. A. Lyman for noting this similarity and pointing 
out that it could be generalized to include the case considered herein. 
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The physical meanings of the various terms in equations (12) are as follows. Since 


[Re(I)] 2 is the time average of the square of the real current density over a cycle, the 
quantity W represents the average rate of joule dissipation; that is, the average power 
expended during a cycle in resistive heating of the plasma. By similar reasoning, Uj^ 
is seen to be the average energy stored in the magnetic field during a cycle. The fore- 
going two terms are those expected for an ordinary solid conductor in a time-varying 
magnetic field. The term , however, is not so intuitively obvious; in fact, it usually 
is not considered in analyses of the r-f heating of solid conductors. It arises in the 
present analysis because of the importance of electron inertia in low-pressure r-f 
plasmas. By writing the real current density as Re( j ) = n e qv e , this term becomes 


U, 


ke 



dV 


which clearly shows that U ke is the time average of the kinetic energy invested in the 
drift motion of the plasma electrons. 2 Both U ke and U M contribute to the plasma 

2 

It is possible to prove in general, that, for a fully ionized plasma with no mean ve- 
locity, U ke is the sum of the time -averaged electron and ion drift kinetic energy. But in 
the present case, it can be demonstrated that the azimuthal ion drift velocity is m e /rm 
times the electron drift velocity, since the plasma as a whole has no motion in the azi- 
muthal direction (note that ( j x B)^ = 0 because j = j = 0). Hence, the ion drift kinetic 
energy is negligible compared with that of the electrons. 
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reactance, as is clear from equation ( 12 a). 

Equation (12a) will be used in the following manner. The complex Poynting vector 
(e x B )/\ Uq will be evaluated at the surface of the plasma, so that the quantity on the 
left hand side of equation (12a) is known. Taking the real part of this quantity will then 
yield the -average resistive power W. The average energy stored in the plasma during 
a cycle U will be obtained from the imaginary part of - J (e x B*/ 2jUo) * dS. An effec- 
tive inductance of the plasma L can then be defined by the relation 

tr 

- LI 2 = U (13) 
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o 

where I is the mean-square current in the r-f coil. 


Solution for Magnetic and Electric Fields 


In order to evaluate the quantity on the left side of equation (12a), it is necessary to 
solve equations (7), ( 8 ), and (9) for B and E. It is convenient to introduce two param- 
eters that will be used throughout this report, namely 
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(15) 


where a is the ratio of the r-f field frequency to the collision frequency v and 5 is 
the normal r-f skin depth. 

The following equation for B is obtained by substituting equations (9) and (14) into 
equation ( 8 ) and taking the curl of both sides 

- i Mn w - 

V x V x B = B (16) 

77 (1 + io) 


Noting that B = Bz and rearranging yield 
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( 17 ) 


d^B + 1 dB i/X()Ct> b - 0 
dr 2 r dr rj(l + ia) 

Using equation (15) and changing variables from r to | where 
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yield Bessel's equation of zero order 


d2 ^ + 1 f® + B = 0 
d| 2 « 


The solution of equation (19) for B finite at | = 0 is 


B - KJ 0 (|) 

Applying the boundary conditions B = B„ at r = a or ij = £ yields 
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To determine E, equation (7) may be rewritten as 


- il (rE) = -iwB 
r dr 


( 21 ) 


( 22 ) 


where E = E0. Substituting equation (21) into equation (22), using equation (18), and 
simplifying yield 
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which may be written as 


0,6 / 2 ^ 4 J 1 <€) 

E = - — (l + or) B — exp 
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(25) 


The variation of E and B with radial position is considered in detail in appendix B . 
The significance of two important limiting cases should be mentioned at this point, how- 
ever. For a = 0, that is, for high density plasmas where the electron inertia term in 
the Ohm’s law is negligible, the solution reduces to that for the induction heating of a 
solid cylindrical rod (ref. 8). Hollister (ref. 9) has investigated this particular case as 
it applies to a diagnostic method for high-density plasmas. In this case, it is well known 
that the inward penetration distance of the electric and magnetic fields is the normal r-f 
skin depth 6. On the other hand, for low -density plasmas a » 1 and from equation 
(18), it follows that, in the limiting case of very large a, 


where 






(26) 


(27) 


In this case, the penetration distance becomes c/u^, which is the characteristic attenu- 
ation distance for a transverse electromagnetic wave in a collisionless plasma when the 
wave frequency is much below the plasma frequency (ref. 6, p. 54). 


Power Dissipated in Plasma 

The solution for the power dissipated in the plasma can now be obtained by evaluating 
the real part of the complex Poynting vector at the surface r = a. When r = a and 
| = £ , the left side of equation (12a) becomes 
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By taking the real part of equation (32), the average power coupled to the plasma per unit 
surface area is 


p.JL 


= ^^^[(l + 0 2 ) 1/2 + »] 1/2 Re 


27ra l 4 p 


Now, define a nondimen sional power P as 
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and the quantity h as 
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Note that h is the magnitude of the complex argument | and is thus the ratio of radius 
to field penetration. 

Combining equations (34) and (35) yields 
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where equation (18) becomes 
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A plot of P against h with a as a parameter is given in figure 1. The cases 
a = 0 (high-density plasma) and a » 1 (low-density plasma) are shown by the curves, 
while some intermediate cases (a = 1, 3) are indicated by data points. Figure 1 shows 
that P is independent of a for a > 3 and that the solution for a » 1 may be used 
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Figure 1. - Nondimensional power flux as function of radius to field penetration 
ratio with a as parameter. 


as a very good approximation for these a values. This solution is shown in appendix C 
to be 


P = 


2V2 


I?(h) I 9 (h) 


In(h) 




( 38 ) 


where the I(h)’s are the modified Bessel functions. In this case, h is a function of 
electron density only; that is, 


h = 



co a 
P 
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( 39 ) 
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Inductance of Plasma and Inductance Change of Radio-Frequency Coil 


The average energy U stored in the plasma during a cycle is obtained by taking the 
imaginary part of equation (32) and dividing by 2o> (see eq. (12a)). Thus, 


- B> 'ai[(i + a 2 ) 1/2 + 0 ] 1/2 im 


u = 
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The effect of this energy storage on the external circuit in an experiment is to reduce 
the observable inductance of the r-f coupling coil. Defining the inductance of the plasma 
such that 



(see ref. 7, p. 202) and noting that 1^ = 1^/2 yield 
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where 

Y = [(l + a 2 ) V2 + a ] 1/2 lm 

Consider the observable change in the coil inductance when the radius of the r-f 
coil b is larger than the plasma radius a. The observed change of inductance will 
be the change that occurs in the volume occupied by the plasma. When no plasma is 
present, the inductance associated with this volume is 
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The change of inductance when a plasma is present is, therefore, 
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AL = L 
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which may be rewritten as 



(45) 


The quantity AL/L (eq. (45)) is plotted as a function of h with a as a parameter in 
figure 2. It is again seen from this figure that the a dependence is not important for 
a > 3 and, therefore, the solution for a » 1 may be used for these a values. In 
this case, the observable change in coil inductance will be a function of electron density 
only and is given by (see appendix C) 


AL 



2 # 

h I Q (h)_ 


(46) 


It should be noted, however, that the fractional change in inductance, which is measured 
experimentally, is AL/L T , where Lrp is the total inductance of the r-f coil. Conse- 
quently, even though AL is a large fraction of L^, the measurement may be quite 
difficult if it is a small portion of L,p. 



Figure 2. * Nondimensional inductance change as function of radius to field 
penetration ratio. 
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EXPERIMENT 


Apparatus 

A photograph of the r-f induction heater in operation is shown in figure 3. Figure 4 
is a schematic diagram of the apparatus. The test chamber consists of two concentric 
glass funnels with cooling air passing between them. Two glass to metal seals are at- 
tached to the top of the inner funnel and have conventional O-ring fittings welded to the 
metal sections. Helium gas was introduced into the system through a quartz tube in- 
serted in the center O-ring fitting. The second fitting was used as an entry port for di- 
agnostic devices. 

A Faraday shield, consisting of a series of insulated vertical copper wires with a 
common ground, was located between the funnels to ensure operation in the inductive 
mode. A small confining direct -current field of 20 to 30 gauss (20xl0 - ^ to 30x10’^ T) is 
supplied by the direct-current field coils in a magnetic mirror geometry. This field in- 
creases the stability of the discharge and allows continuous operation for long periods of 
time. 

The 7. 7 -centimeter -diameter r-f coil is a 4-turn silver-plated copper coil of length 
7. 5 centimeter. Magnetic -field measurements were made with a single turn magnetic- 
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Figure 4. - Schematic of test apparatus. 


field probe. This probe was located on the wall of the inner funnel; hence, it measured 
an effective average of the r-f field over the cross sectional area. Measurements at 
various axial positions showed that the r-f magnetic field was approximately constant 
(±10 percent) over the coil length. Both the direct-current and r-f coils were water 
cooled. The r-f power is supplied by a 2 -kilowatt, 17. 5-megahertz transmitter. The 
capacitors and Cg are calibrated, variable vacuum capacitors. 

A microwave interferometer is used to measure electron density. Two microwave 
windows were positioned on the sides of the inner funnel. A fiber optics light pipe was 
used to transmit the helium spectral line radiation to a 0. 5-meter Ebert monochromator. 
The entire system was exhausted into a 32-inch (0. 813 m) oil diffusion pump. 


Electrical Instrumentation and Determination of Plasma Properties 

The electrical properties needed to compare the experimental data with the theory 
are the power coupled to the plasma, the current in the r-f coil, the change of inductance 
of the r-f coil, and the plasma electron density and resistivity. 
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Power measurement. - The power coupled to the plasma is measured by using a 
commercial directional coupling wattmeter. This wattmeter measures the forward power 
into the matching network and load coil when the system is tuned for zero reflected power 
and its impedance matches the transmitter impedance of 50 ohms. 

Coil current and inductance measurements. - The values of the coil current and in- 
ductance change are obtained from direct measurements of the circuit parameters. The 
equivalent circuit for the impedance matching network, coil, and plasma system is 
shown in figure 5. 



In this figure, Eq is the applied voltage and Rq is the total impedance when the 
system is tuned for zero reflected power and equals the transmitter impedance of 
50 ohms. C l and C 2 are calibrated variable capacitors and L 2 is the inductance of 
the strap connecting Cj and C 2 plus the lead inductance of the coil. L is the coil induc- 
tance plus the lead to ground inductance, R is the coil resistance, and C g is the stray 
capacitance. This capacitance is a result of the distributed capacitance of the coil, the 
Faraday shield, and the voltmeter. However, because the effect of this capacitance on 
the circuit measurements is the main interest, it is considered as a lumped capacitor 
across the r-f coil. The plasma inductance and resistance are represented by Lp and 
R , respectively. The ammeters A 2 and Ag measure the rms values of the current 
I 2 and the coil current Ig, while the voltmeter V p measures the peak voltage across 

the L-R-C network. For a high Q circuit, the change in inductance may be obtained 
s 

directly from the V and Ig measurements and the relation 
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( 47 ) 


L - AL = — — - 

V2 I 3 o, 


However, because of some inconsistencies in the voltage measurements (see appendix D), 
an alternate method for determining AL, which did not require voltage measurements, 
was also employed. The net reactance of the L - AL, R + AR, C portion of the circuit 
is inductive, and this portion may be replaced by an apparent resistance R„ and appar- 

a 

ent inductance L a ^. For a high Q circuit, the magnitude of the impedance of this por- 
tion of the circuit is, to a very good approximation, that of the apparent inductance (see 
Swett, ref. 10); that is, 


J eff 


(r 2 + » 2 l 2 ) 1/2 


w 2 C 2 R 2 + (o> 2 LC g - l) 


n 1 / 2 



(48) 


Therefore, the voltage V is also given by 


V p - Vi I 2 u,L ap 


and using equation (47) yields for L - AL 


L " AL = L ap 
x 3 


(49) 


(50) 


At resonance conditions (Zq = Rq = 50 ohms), L a p may be obtained from the and Cg 
readings using the relation 


P R 2 

_ C 1 K 0 , 1 _ T 
a P 2 2 2 2 2 

w C^Rq + 1 w z C 2 


or using the values Rq = 50 ohms and o> = 1. 099*10° per second results in 
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(51) 


(52) 
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where and Lg are in microhenries if Cj and Cg are given in picofarads. 

In order to ensure the validity of the current and voltage measurements and the 
values of the circuit parameters, it was necessary to conduct a detailed circuit check 
(presented in appendix D). 

Plasma property measurements. - The electron density was obtained by measuring 
the phase change of a microwave beam when it passed through the plasma. This mea- 
surement was accomplished by using the 39. 5 gigahertz microwave interferometer that 
is shown schematically in figure 6. The data was reduced following the method of Kuhns 
(ref. 11). 
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Figure 6. - Schematic diagram of microwave interferometer. 


The plasma resistivity is obtained from measurements of the pressure, electron 
temperature, and electron density and by using the relations given in appendix E. The 
electron temperature is obtained by using the spectroscopic technique refined in refer- 
ence 12, which is also discussed in more detail in appendix B. 


Experimental Procedure 

The procedure followed was to first obtain data for the system with no plasma pres- 
ent. In these tests, the impedance matching network was tuned for zero reflected power 
and the forward power was measured for a series of coil currents. Once these measure- 
ments were obtained, the desired helium gas flow was set and the r-f voltage was then 
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increased until a discharge was initiated. Once the plasma was obtained, the match net- 
work was again tuned to zero reflected power and the desired coil current was obtained 
by varying the r-f voltage. The system was then allowed to run for a few minutes, and, 
if no changes in the reflected power or coil current were observed, the data were taken. 

The quantities measured for each data point were the pressure, forward power, coil 
current, direct field current, transmitter -plate voltage and current, series and shunt 
capacitances, spectral line intensities, microwave interferometer measurements, and 
the various circuit voltages and currents needed to determine the inductance change and 
other circuit parameters. 


Experimental Data 


Some typical experimental data are presented in table I. In obtaining the experi- 
mental data, the power coupled to the plasma P , the pressure, electron density, coil 
current, and electron temperature were found to be interdependent quantities. Because 
of the limitations on the amount of power (0 to 1000 W) that could be coupled to the 
plasma and the range of electron densities (1x10 to 5x10 m” ) that could be mea- 
sured, a detailed parametric study of these quantities was not made. However, a typical 
variation of these quantities with pressure at a constant coil current is shown in figure 7. 
This figure shows that the electron density and the power coupled to the plasma increase 
almost linearly with pressure. The electron temperature decreased with increasing 
pressure, a result which is in agreement with the results of Pell for low pressures 
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Figure 7. - Typical variation of power coupled 
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TABLE I. - DATA TABLE 


Root mean Forward 
square coil power, 
current, W 

A i 

Total power 
coupled to 
plasma, 

*g- 

W 

Power transfer 
efficiency, 
percent 

Electron kinetic 
temperature, 
kT e> 

eV 

Electron 

density, 

V 

m- 3 

Neutral 

density, 

V 

m' 3 

Fraction 

ionized, 

f. 

i 

Plasma 

resistivity, 

V, 

ohm-m 

m' 1 

Inductance 

change, 

AL, 1 

Method 1 

Method 2 

14 

219 

108 

50 

8.8 

6. 7X10 16 

5. 34xl0 20 

1. 25xl0' 4 

19 xlO' 3 

14. 5 

0.03 


14.8 

247 

127 

51 

8.7 

7. 1 

5.05 

1.41 

17 

17 

.039 

0. 043 

15.6 

289 

152 

53 

9. 1 

10. 2 

5.0 

2. 04 

11 

27.9 

.031 

.036 

16.4 

335 

182 

55 

8.4 

11.2 

4. 97 

2.25 

10 

33.2 

.043 

.027 

16.9 

490 

325 

66 

4.0 

36.7 

7.39 

4.96 

4.8 

119 

.045 

.068 

17.4 

423 

253 

60 

7.3 

17.1 

5.42 

3. 15 

7.0 

58 


.018 

18 

390 

195 

50 

— 

12 

— 

— 

--- 


.029 

.014 

18 

364 

169 

47 

— 

6.5 

____ 

— 

--- 


.018 

.oi ; 

18 

474 

279 

59 

— 

19. 5 

— 

— 

--- 


.036 

.033 

18 

510 

315 

62 

— 

22.8 

— 

— 

— - 


.036 

.040 

18 1 

575 

380 

66 

— 

35. 5 

— 

— 

— 


.051 

.051 

18 

620 

423 

68 

1 

40 

— 

— 

--- 


.047 

.067 

18. 1 

435 

240 

55 

8. 5 

: 14.9 

5.0 

2.98 

8 

48 

.036 

.023 

18. 1 

594 

404 

68 

4.0 

1 42.8 

7.4 

5.8 

4.3 

144 

.056 

. 070 

18.3 

474 

279 

59 

6.8 

; 19.9 

5. 34 

3.72 

6. 6 

63 

.030 

.025 

19.2 

505 

300 

60 

8.3 

j 17. 1 

4.79 

3. 56 

6. 5 

63 

.045 

.040 

19.7 

705 

480 

68 

4.0 

j 49 

7.08 

6.92 

3.4 

182 

.066 

.068 

20 

l 

548 

323 

59 

8. 1 

! 18.9 

4.73 

4.0 

6.0 

71 

.044 

.052 ! 

20 

594 

372 

63 

6.5 

124 

5.4 

4.45 

5.4 

86 

.045 

.049 ! 

20 

621 

396 

64 

6. 1 

: 26.5 

5. 58 

4.75 

5.2 

95 

.052 

.048 

20 

715 

475 

67 

— 

' 33. 2 

— 


— 


.036 

. 064 

20 

760 

520 

69 

--- 

44 

— 

— 

--- 


.046 

.068 

20 

895 

655 

73 

! 

; 42. 5 

— 

— 

--- 


.042 

.062 ; 

| 20. 1 

588 

363 

62 

1 6.7 

: 23 

5. 23 

4.4 

5.4 

: 83 

.043 

. 036 ! 

! 20. 5 j 

i 

622 

382 

62 

6.7 

19. 5 

5. 26 

3.7 

6.6 

69 

.047 

. 070 

i 

20.6 

676 

435 

64 

5.9 

24. 1 

5.65 

4. 26 

5.6 

92 

.052 

.062 

21. 2 j 

715 

452 

63 

6.0 

28. 1 

5.39 

5.21 

4.5 

; 112 

.045 

.047 

21.8 

695 

425 

61 

7.6 

23.8 

4. 66 

5. 10 

4.3 

' 104 

.043 

.042 : 

22. 1 

798 

523 

j 66 

j 6.0 

' 32 

5. 34 

6.0 

4.0 

| 134 

.040 

.063 




(private communication from K. M. Pell of University of Florida). The power transfer 
efficiencies defined as the power coupled to the plasma divided by the forward power are 
also listed in table I. These efficiencies were from 47 to 73 percent. 


COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 

The solution for the nondimen sional power per unit area coupled to the plasma for 
the case a » 1 was given by equation (38) as 

Ij(h) I 2 (h) 

I 2 (h) I 0 (h) 



where h = co^a/c = 1.88x10”^ therefore, P is a function of n 0 only. 

The dimensional power coupled to the plasma per unit area is obtained from equa- 
tion (34) and is 



Noting that 


sin 



a 


— for a » 1 yields for P 
2a 


_ B 2 wa _ 

P = — 1 P(n e ) (54) 

V2 J UqQ! 

The total power coupled to the plasma is, therefore, 

P g = 27raZP = 27ra2 ^ B 2 P(n e ) (55) 

V2 m 0 « 

In order to compare the results of this analysis with the experimental data, the following 
substitutions are made: 
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and 


a = 


m a u 

e 


r?n e q 


n 2 , 2 , t 2 t 2 
2 C b m 0 N I 

B a = r- 


m 


where C B is a correction factor for the short solenoid used in the experiment. The 
value of C B was obtained from measurements of the field produced by the coil at differ- 
ent axial positions with magnetic -fie Id probes. The following values of the experimental 
constants were used: 


Plasma radius, a, m 

Length of plasma, l, m 

2 

Correction factor for short solenoid, Ctj 

-1 D 
Applied frequency, u>, sec 

Number of turns in coil, N 


. . 0.024 
. . 0.075 
. . 0.233 
1.099X10 8 
4 


Substitution of these values and the quantities 


a and B 


2 

a 


into equation (55) yields 


P g = 4.52xi0" 15 4r7n e P(n e ) 


(56) 


or 



4. 52X10' 15 n P(n ) 

v C 


(57) 


The assumption a » 1 produces a maximum error in P of 4 percent if a = 3 
and 2 percent if a = 4. Most of the a values obtained experimentally are a > 4; 
therefore, the experimental results are compared with the theory for a » 1. In this 
case, AL is given by (see appendix C) 


AL = 



2 

h yhj_ 


( 46 ) 


where 



is defined by equation (43) and is calculated to be 0. 113 microhenry. 
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Figure 8. - Comparison of experimental and theoretical 
results for power coupled to plasma (see eq. (57)) for 
a » 1. 



Figure 9. - Comparison of experimental and theoretical 
results of inductance change for a» 1. 


A comparison of the theoretical (eqs. (57) and (46)) and experimental results is 
shown in figures 8 and 9. Reasonably good agreement is evident between theory and ex- 
periment. The experimental data for the power coupled to the plasma are slightly higher 
(«10 percent) than the theory predicts, but they definitely exhibit the same dependence on 
electron density that is predicted by the theory. One reasonable explanation for the 
slight discrepancy in magnitude could be an error in the determination of the plasma 
volume. The plasma radius was taken to be the inner radius of the inner funnel, and the 
plasma length was taken to be the length of the r-f coil. In reality, the plasma did fill 
the coil and inner funnel but there was also a plume of plasma above the r-f coil. The 
tendency would be to underestimate slightly the plasma volume. This difference could 
easily account for the slight discrepancy between the theoretical and experimental re- 
sults. Another possible reason for this discrepancy could be in the treatment of the 
electron density. In the theory, the spatially dependent electron density was replaced 
by a uniform effective electron density. The microwave interferometer measures an 
average electron density across the diameter of the plasma. However, this averaged 
electron density is not necessarily the same as the assumed effective electron density. 

The direct-current field did not affect the data since the system was not operated 
at any of the resonant frequencies for the plasma constituents. Data taken with and 
without the direct -current field fell on the same curve. However, the direct -current 
field permitted steady-state operation over a wider range of parameters. 

The inductance change data (fig. 9) also agreed with the theory, although there is a 
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large amount of scatter in the data. This scatter is understandable, however, in view of 
the difficulties entailed in measuring such small inductance changes. It was for this rea- 
son that two methods of measuring these changes were employed. An experiment de- 
signed such that a large change in inductance is expected would probably have much less 
scatter. 


CONCLUDING REMARKS 

An experimental and theoretical investigation of radio -frequency induction heating 
and production of plasmas has been carried out. A linearized analysis was performed 
for the coupling of a solenoidal radio -frequency magnetic field with a steady cylindrical 
plasma. The good agreement obtained between the experimental and theoretical results 
indicates that the linearized analysis serves as a very good approximation for predicting 
the power coupled to the plasma and inductance change as a function of the plasma prop- 
erties. The fact that the inductance change is a function of electron density only for low- 
pressure plasmas indicates that the measured inductance change may be used to measure 
the average plasma electron density. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 25, 1967, 

120-26-03-04-22. 
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APPENDIX A 


SYMBOLS 


A 2’ A 3 
A' 


a 

B 

B’ 



c v c 2 

c 


D 

d’ 

E 

E 0 

f i 

H 

h 

I 

Im 


*n< h > 

h 

X 3 

J n (?) 


ammeters 

2. 3 log 10 d'/c 

plasma radius, m 
magnetic -flux intensity, T 
0.366 (1 - 2.5/N+ 3.8/N 2 ) 
coil radius, m 
capacitance, pF 

correction factor for magnetic field produced by short coil 

correction factor used in calculating coil inductance, defined by eq. (C3) 

stray capacitance, pF 

calibrated variable capacitators 

winding pitch of coil 

2 

electric displacement, c/m 
diameter of coil tubing, m 
electric -field strength, V/m 
applied voltage, V 
fraction ionized 
magnetic -field intensity, A/m 

nondimensional radius to field penetration ratio (see eq. (35)) 
electric current, A 
imaginary part of 

modified Bessel function of order n and argument h 
current through Cg, A 
coil current, A 

Bessel function of order n and argument £ 



j ratio of rms magnitude of magnetic field at some radial position r to that at 

boundary r = a (eq. (B4)) 

ratio of rms magnitude of electric field at some radial position r to that at 
boundary r = a (eq. (B12)) 
o 

j current density, A/m 

K constant defined by eq. (20) 

K n correction factor for short coils defined by eq. (C2) 

k Boltzmann constant, 1. 3805xl0 - ^ J/°K 

kT g electron kinetic temperature, eV 

L inductance, pH 

Lq coil inductance, pH 

L Q p inductance of plasma volume when no plasma is present, pH 
Lp plasma inductance, pH 

L 2 inductance of strap used to connect Cj and C 2 plus lead inductance of coil, pH 

l length of plasma, m 

M n (h) modulus of Kelvin functions of order n and argument h 

Q 1 

m e electron rest mass, 9. 108x10 kg 
N number of turns in coil 

_Q 

n g electron density, m 

_ q 

nQ neutral density, m 

P average flux to plasma during cycle, W/(27ra l), W/m^ 

P total power coupled to plasma, W 

o 

P nondimensional power defined by eq. (34) 

p* pressure, N/m or pm 

Q coL/R 

-IQ 

q electron charge, 1.602x10 C 

R resistance, ohm 

R apparent resistance, ohm 

Re real part of 
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R q total impedance when system is tuned for zero reflected power, 50 ohms 

Rp plasma resistance, ohm 

r radial coordinate 

2 

5 surface area, m 

T temperature, °K 

Tq neutral temperature in test chamber 

Tq neutral temperature in McLeod gage 

t time, sec 

U average reactive energy in plasma during cycle 

U^ e average drift kinetic energy of electrons during cycle 

U M average magnetic energy during cycle 

3 

V volume, m 

V peak voltage, V 

r 

v velocity, m/sec 

v e electron velocity, m/sec 

W average resistive power coupled to plasma during cycle (see eq. (12b)) 

X reactance, ohm 

V defined by eq. (39) 

Z impedance 

Z eff effective impedance of L-R-C g portion of circuit 
z axial coordinate 

a dimensionless coordinate defined as ratio of r-f field frequency to collision 

frequency 

/3 tan" 1 [rhl 1 (rh/a)/(2aal 0 (rh/a)j, see eq. (B8) 

0' tan _1 [hl 1 (h)/2o'l 0 (h)], see eq. (B8) 

y tan” 1 [2ao'I 1 (rh/a)/rhI 2 (rh/a)], see eq. (B16) 

y' tan -1 [2al 1 (h)/hl2(h)], see eq. (B16) 

6 normal r-f skin depth, (2rj/ 
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e 0 

V 

^ea 

^ei 

6 

e n (h) 

A 


M 0 



ea 


i 

P 


a 


m 


el 


dielectric constant of free space, F/m 
plasma resistivity, ohm-m 

resistivity due to electron -atom collisions, ohm-m 
resistivity due to electron-ion collisions, ohm-m 
azimuthal coordinate 

phase angle of Kelvin functions of order n and argument h 

12 7r(e 0 kT e ) 3/2 /q 3 ny 2 

magnetic permeability of free space, H/m 

total electron collision frequency, sec”* 

electron-atom momentum transfer collision frequency, sec”* 
v e2L / n 0 > m 3 /sec 

dimensionless variable defined as (r/6)[-2i/(l + io: )]l/2 
ratio of radius to r-f skin depth, a/6 

o 

total momentum transfer cross section due to elastic collisions, m 


t time between electron collisions or total electron collision time, sec 

cp tan”* Xg/R a 

Sf electron kinetic energy, eV 

co applied frequency, rad/sec 

a> r resonant frequency of L-C network, rad/sec 

Subscripts: 

a evaluated at radial position r = a 

ap apparent 

m maximum value 

Superscripts: 

~ time averaged value over cycle 

* complex conjugate 
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APPENDIX B 


RADIAL PENETRATION OF MAGNETIC-FLUX DENSITY AND 
ELECTRIC-FIELD STRENGTH 


In this appendix, the variations of B and E with radial position are investigated, 
and tables and plots of these variations are presented. 


Solutions for Magnetic-Flux Intensity 

General solution . - Including the time dependence, the solution for B is given by 
(see eq. (21)) 


B = B Re 

ct 


J 0^> e iwt 

lW 6 - 


(Bl) 


where 


and 


V 2 


| = — exp 

5 / o\l/4 


(l + a 2 ) 


Al + icot^a\ = r h fix + i co CM 
\2 2 ) a \2 2 / 


t a -i exp, 

a 6 , o\l/4 \2 


(l + a 2 ) 


(in + i cof 1 ^ = h e /irr + i cot/V 

2 / V 2 2 J 


The quantity J 0 (|)/J 0 (i; a ) is a complex number which may be written in the polar form 

J 0<» _ 


J 0<«a> 


) 


(B2) 


Equation (Bl) then becomes 


B = B a ^Re £ e i(wt+<p)j = C os(a)t + cp) 


(B3) 
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Taking the root mean square value yields 


B 


rms 


( B a), 


rms 


(B4) 


Therefore, the ratio of the root mean square magnitude of B at some radial position r 
to that at the boundary r = a is given by j , and the phase difference between the two is 
given by cp. 

Solutions wh en a = 0 and a » 1. - When a = 0, 

{ = £ hi 3 / 2 
a 


and 


k - hl3/2 


yo _ M °(; h ) 

ty«a> M 0 (h > 



- i^ 0 (h) 


(B5) 


where M Q and are the modulus and phase, respectively, of the Kelvin functions. 
Therefore, when a = 0, 


and 



M °(a h ) 
M 0 (h) 



- 


When a » 1, 

{ = £ hie 1/2 “ 
a 


and 


= hie i/2a 
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It is also shown in appendix C that, for a » 1, 


J o<« ' l o(l h ) 



and 


J 0 (5 a ) = IqW + A U!<h) 


Therefore, Jg(£)/JQ(£ a ) in polar form is, for a » 1, 

.2 o \”il/2 




J 0<«a> 


!#»♦(£)>] 


1/2 


e i(/3-/3’) = j e i(P-(3') 


(B6) 


(B7) 


(B8) 


where 


and 


/3 = tan 



/3 1 = tan 


hlj(h) “ 
2o!l 0 (h) 


For a » 1, and from equation (B8), 



and 

<P = P~I 3' 
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Solutions for Electric-Field Strength 


General solution . - Evaluating equation (25) at r = a yields, for E including the 
time dependence, 


E = E„Re 

d 


J x (l) 



(B9) 


Following the procedure used in the preceding discussion for the magnetic field results 
in 


J x (i) 



and 


E = E a /i cos(cot + (pj 


(BIO) 


(Bll) 


or 


rms 


( E a) 


= h 


rms 


(B12) 


Therefore, the ratio of the rms magnitude of E at some radial position r to that at the 
boundary r = a is given by f and the phase difference between the two is cpy 

Solutions when a = 0 and a » 1 . - Letting a - 0 and using the results of the 
first section of this appendix reveal that 


and 


_M^ 

Mj(h) 



When a » 1, it is shown in appendix C that 


(B13a) 


(B13b) 
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r 


and 


jj<«) = el/2 “[ u i(f h ) 



J,(«J = e i/ 2 “[il 1 (h) - A I -(h) 
1 a L 1 2a * J 


Therefore, in polar form, 


jjd) 


P(;‘H 

rh > 
2 a a, 

2 0 
) i\ 

P 14 

■j 1/2 


f h ) 

\2 a) 

2 o 

i ■!<>•> 

1/2 


e i(y-r') 


where 


and 



r' 


= tan 


” hI 2 ( h ) _ 


For a » 1, equation (B16) yields 


h- 



<P 1 = y - y' 


(B14) 


(B 15) 


(B16) 


(B17a) 

(B17b) 
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Tabulation and Use of Results for Radial Variation of Magnetic-Flux 
Density and Electric-Field Strength 


Numerical values of the rms magnetic- and electric -field ratios } and re- 

spectively, are presented in table n as functions of r/a, h, and a. The values for the 
phase angles cp and cp^ are given in table HI. In table HE, it is understood that the 
fields in the plasma lag behind those at the boundary r = a and the absolute values of the 
phase angles are given. 

In table II, results are presented for a = 0, 1, 2, and a » 1 only because the 
a » 1 solution is good for a > 3. When using results of this type the reader is gener- 
ally interested in plots of the radial variation of the fields ( J, for a given r-f skin 

depth and various a values. These results may be obtained from the data presented in 
table n. The manner in which the data are presented in table n is advantageous in that 
the data encompass a wide range of 5 and a values. Letting the ratio of the radius to 
the r-f skin depth be p (i. e. , p = a/6) yields for h 


h = 


(i 


Vi p 
2V/ 4 

f a J 


(B18) 


K / and ^ are plotted as a function of h with r/a as a parameter for fixed a 
values of 0, 1, 2, and a » 1 the desired information may be obtained directly from 
these plots. One merely chooses the desired p value and calculates the h for each a 
needed. The variations of j and ^ with r/a are obtained from the family of curves 
for the chosen a. An example of the type of results obtained is shown in figure 10 for 


p = 10. 

The data for the phase angles cp and cp^ are presented in table HI. The phase 
angles depend on a even when a » 1; therefore, these results are presented for a 
values of 0, 1, 2, 3, 5, 25, and 100. In order to obtain data for a constant p value, the 
same technique used for the J and j ^ quantities is followed (i. e. , a family of curves 
of cp and cp ^ as a function h is obtained with r/a as a parameter for each a). If 
data for a values not included in table III are needed, the data in table in must be used 
to obtain auxiliary plots of cp and cp^ as a function of a with h as a parameter for 
each r/a. A plot of the phase-angle variations with r/a and a for p = 10 is pre- 
sented in figure 11. 
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TABLE n. - VARIATION OF ROOT MEAN SQUARE MAGNITUDES OF MAGNETIC-FLUX AND 
ELECTRIC- FIELD STRENGTH WITH RADIAL POSITION 


[Radial position r normalized to 1 at boundary r = a!] 


Nondi mensional 

Root mean 



Nondi mensional radius to field penetration ratio, 

h 


radius, 

square 

- - 

i 


\ 

r 

1 

| 

■ 

— 



r/a 

magnitudes 

0. 1 

0.5 

1.0 

2.0 

2.5 

| 5.0 

1 7 * 5 

[ 10 

15 

20 

25 






a 

= 0 







0 


0.999 

0.998 

0.983 

0. 812 

0. 662 

0. 160 

0. 034 

0. 007 

2xl0' 4 

8xl0" 6 

5X10’ 8 


h 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.1 

i 

0.999 

0. 999 

0.984 

0. 814 

0. 66 

0. 161 

0.034 

0. 007 

3xl0“ 4 

1X10“ 5 

-7 f 

4X10 ' 


i\ 

. 100 

. 100 

. 100 

.096 

.091 

.043 

.013 

.003 

2X10" 4 

9xl0' 6 

4xl0" 7 i 

0.2 


0. 999 

0. 999 

0. 984 

0. 814 

0. 662 

0. 163 

0. 036 

0. 008 

5xl0 r4 

3x10’ 5 

2xl0“ 6 

l 

h 

. 200 

. 200 

. 199 

. 192 

. 182 

.086 

.027 

. 007 

4xl0' 4 

3xl0’ 5 

2x10’ 6 j 

] 0. 3 

f 

0.999 

0.999 

0.984 

0. 815 

0. 664 

0. 173 

0. 046 

0. 013 

lxlO' 3 

9xl0’ 5 

-6 1 
8X10 D 

i - 

/ 1 

.300 

. 300 

. 299 

. 288 

. 273 

. 131 

, 042 

.012 

lxlO' 3 

9xl0’ 5 

8X10' 6 

0.4 

/ 

0.999 

0. 999 

0.984 

0.819 

0. 672 

0. 197 

0,066 

0. 023 

3x10’ 3 

3xl0~ 4 

4x10' 5 


h 

.400 

.400 

. 399 

. 385 

. 365 

. 179 

.064 

.022 

3x10’ 3 

3xl0" 4 

4x10' 5 

0. 5 

i 

0.999 

0. 999 

0.985 

0. 827 

0. 686 

0. 243 

0. 100 

0. 042 

7xl0~ 3 

0.001 

2xl0" 4 


f i 

. 500 

. 500 

.499 

. 482 

.458 

. 236 

.097 

.040 

-J 

X 

o 

1 

.001 

2xl0' 4 

0.6 

/ 

0. 999 

0. 999 

0.986 

0. 840 

0.712 

0. 313 

0. 156 

0. 077 

-2 

1.9x10 z 

0.005 

lxlO' 3 


K 

. 600 

.600 

. 599 

. 579 

. 553 

. 310 

. 151 

. 075 

1. 8x10' 3 

.004 

lxlO' 3 

0.7 

/ 

0. 999 

0. 999 

0.988 

0, 862 

0. 753 

0.414 

0. 244 

0. 143 

5X10“ 2 

0.017 

4xl0" 3 


h 

.700 

.700 

. 699 

. 679 

. 652 

.419 

. 240 

. 142 

4.9X10’ 2 

.017 

6xl0' 3 

0.8 

i 

0.999 

0.999 

0.990 

0. 893 

0. 813 

0. 55 

0. 388 

0. 272 

1. 34x10" 1 

0.066 

3. 3x10" 2 



.800 

.800 

. 799 

. 782 

. 758 

. 546 

. 384 

. 270 

1. 33x10“ 1 

.066 

3. 2x10’ 2 

0.9 

} 

0.999 

0. 999 

0.994 

0.939 

0. 895 

0.74 

0. 62 

0. 520 

3. 65x10" 1 

0. 256 

1. 8x10' 1 



.900 

.900 

. 899 

. 888 

. 873 

. 736 

. 618 

. 518 

3. 64x10“ 1 

. 260 

1. 8x10" 1 

1.0 

t 

1.000 

1. 000 

1. 000 

1.000 

1.000 

1.000 

1. 000 

1.000 

1.000 

1 

1.000 

1.000 



1.000 

1.000 

1.000 

1.000 

1.000 

1. 000 

. i 

1.000 

1.000 

1. 000 

1.000 

1.000 
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TABLE n. - Continued. VARIATION OF ROOT MEAN SQUARE MAGNITUDES OF MAGNETIC -FLUX AND 
ELECTRIC- FIELD STRENGTH WITH RADIAL POSITION 


[Radial position r normalized to 1 at boundary r = a.J 


Nondimensional 

radius, 

r/a 

Root mean 



Nondimensional radius to field penetration ratio, 

h 


square 

magnitudes 

0. 1 

[ 0.5 

,0 

J 2.0 

| 2.5 

5.0 

7.5 

10 

15 

| 20 

25 






a 

= 1 








0 

/ 

0. 998 

0.956 

0. 836 

0. 510 

0. 370 

0.054 

0.007 

8xl0" 4 

9x10® 

lxlO" 7 

lxlO" 9 


/l 

0 

,0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 1 

/ 

0. 998 

0.957 

0.840 

0. 517 

0.375 

0.056 

0. 007 

9xl0" 4 

lxlO" 5 

2xl0" 7 

3xl0" 9 


A 

. 100 

.098 

.092 

.071 

.059 

.015 

.003 

4xl0" 4 

9x10® 

2xl0" 7 

3x10" 9 

0. 2 

/ 

0.998 

0.958 

0. 845 

0. 528 

0. 387 

0. 064 

0.010 

0.001 

4x10" 5 

9X10' 7 

2x10" 8 


A 

* 200 

. 196 

. 184 

. 143 

. 119 

.033 

.006 

.001 

3xl0" 5 

8xl0' 7 

2xl0' 8 

0. 3 

/ 

0. 998 

0.961 

0.852 

0. 547 

0.409 

0. 079 

0.015 

0. 003 

lxlO -4 

4x10"® 

2xl0" 7 


A 

.300 

.294 

. 277 

. 218 

. 184 

.054 

.012 

.003 

lxlO" 4 

4x10"® 

to 

X 

t— *■ 

o 

1 

—a 

0.4 

/ 

0. 999 

0.964 

0.863 

0. 575 

0.442 

0. 105 

0.025 

0. 006 

4xl0" 4 

2x10' 5 

2x10" 6 


A 

.400 

.393 

.371 

. 298 

. 255 

.084 

.023 

.006 

4xl0” 4 

2X10" 5 

1x10" 6 

0. 5 

/ 

0.999 

0.967 

0. 877 

0.612 

0.487 

0. 145 

0.045 

0.014 

0. 001 

lxlO -4 

lxlO" 5 


A 

.500 

.492 

.468 

. 385 

. 334 

. 128 

.042 

.013 

.001 

1X10“ 4 

lxlO" 5 

0.6 

/ 

0. 999 

0.972 

0. 894 

0.660 

0. 547 

0. 208 

0.082 

0.032 

0.005 

8xl0" 4 

lxlO" 4 


A 

.600 

. 592 

. 567 

. 480 

.426 

. 192 

.078 

.031 

.005 

8x10 - 4 

lxlO' 4 

0.7 

> 

0. 999 



0. 978 

0. 915 

0.721 

0. 624 

0. 303 

0. 151 

0. 075 

0.019 

0. 005 

0. 001 


h 

.700 

.692 

.669 

. 586 

. 534 

.288 

. 146 

.074 

.019 

.005 

.001 

0.8 

f 

0.999 

0.984 

0.939 

0. 796 

0.722 

0. 477 

0. 281 

0. 177 

0. 070 

0.028 

0.011 


A 

.800 

.794 

.775 

. 706 

. 661 

.435 

. 276 

. 175 

.070 

.028 

.011 

0.9 

/ 

0.999 

0.992 

0.967 

0. 889 

0.846 

0.666 

. 

0. 528 

0.419 

0. 264 

0. 166 

0. 105 


A _ 

.900 

.896 

.885 

. 843 

.815 

. 658 

. 524 

.417 

.263 

. 166 

. 104 

1.0 ; 

/ 

1.000 

1.000 

1.000 

1. 000 j 

1.000 

1.000 

1.000 

1.000 

1. 000 

1. 000 i 

1. 000 

i 

L 

/i 

1.000 ! 1.000 

1. 000 

1. 000 j 

1. 000 | 

1.000 

1.000 

1.000 

1.000 

1. ooo 1 

i 

1. 000 
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TABLE n. - Continued. VARIATION OF ROOT MEAN SQUARE MAGNITUDES OF MAGNETIC-FLUX AND 
ELECTRIC- FIELD STRENGTH WITH RADIAL POSITION . 


[Radial position r normalized to 1 at boundary r = a[] 


Nondimensional 

radius, 

r/a 

Root mean 
square 
magnitudes 



Nondimensional radius to field penetration ratio, 

h 


0. 1 

0. 5 

1.0 

2.0 

2. 5 

5.0 

7. 5 

10 

15 

20 

25 






a 

= 2 







0 

/ 

“I 

0.998 

0. 945 

0. 807 

0.466 

0. 326 

0. 042 

0.005 

5xl0' 4 

4x10® 

4xl0~ 8 

4x10" 10 


A 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 1 

/ 

0.998 

0.947 

0. 809 

0. 468 

0.330 

0. 044 

0. 005 

6x10" 4 

7xl0” 6 

9X10' 8 

lxlO' 9 


A 

.099 

.097 

.090 

.066 

.053 

.012 

.002 

3xl0" 4 

4xl0" 6 

6x10' 8 

8x10' 10 

0. 2 

/ 

0. 998 

0,-948 

0.814 

0.480 

0.344 

0.052 

0. 007 

0. 001 

2x10'® 

4xl0' 7 

8x10" 8 


A 

. 200 

. 195 

. 180 

. 133 

. 108 

.026 

.005 

7xl0" 4 

2x10' 5 

4xl0' 7 

7xl0' 9 

0.3 

/ 

0.998 

0.951 

0.823 

0. 502 

0. 368 

0. 067 

0.012 

0.002 

7X10"® 

2xl0' 6 

7xl0" 8 


A 

. 300 

.292 

. 271 

. 204 

. 168 

.045 

.009 

.002 

6x10' 5 

2x10' 6 

7xl0" 8 

0.4 

/ 

0.998 

0.955 

0. 836 

0. 533 

0.404 

0.091 

0.020 

0.005 

3xl0" 4 

1x10® 

7xl0' 7 


A 

.400 

. 391 

.365 

. 281 

. 235 

.072 

.018 

.004 

2xl0" 4 

lxlO" 5 

7X10' 7 

0.5 

/ 

0. 998 

0.959 

0.853 

0. 574 

0.452 

0. 129 

0. 038 

0.011 

0. 001 

8x10" 5 

7x10' 6 


A 

. 500 

.490 

.460 

.365 

.312 

. Ill 

.035 

.010 

9X10" 4 

8xl0" 5 

7x10" 6 

0.6 

/ 

0.999 

0.965 

0. 873 

0. 628 

0. 516 

0. 189 

0.071 

0.027 

0.004 

5xl0~ 4 

8x10'® 


A 

. 600 

j .589 

. 559 

.460 

.402 

. 172 

.067 

.026 

.004 

5xl0" 4 

8x10"® 

0.7 

/ 

0. 999 

0. 972 

0. 898 

0. 695 

0. 598 

0. 282 

0. 135 

0.065 

0.015 

0.003 

8xl0' 4 


A 

. 700 

.690 

. 662 

. 567 

. 511 

. 266 

. 131 

.063 

.015 

.003 

8xl0" 4 

0.8 

/ 

0.999 10.980 

0.928 

0. 778 

0.703 

0. 426 

0. 261 

0. 160 

0.060 

0.023 

0. 009 


A 

.800 

.792 

.769 

.690 

.642 

.413 

. 256 

. 158 

.060 

.022 

.009 

0.9 

/ 

0.999 

0.990 1 0.961 

0. 878 

0.835 

0. 650 

0. 509 

0. 399 

0.245 

0. 151 

0. 093 


A 

.900 

.895 

.882 

! .833 

.803 | 

.641 

. 505 

. 397 

. 244 

. 150 

.092 

1.0 

/ 

1.000 

1.000 

1. 000 

1.000 

1. 000 

1. 000 

1. 000 

1. 000 

1.000 

1. 000 

1.000 


/i 

o 

O 

O 

1. 000 

1. 000 

1.000 

1. 000 

1. 000 ! 1. 000 

1.000 

1. 000 

1.000 

1.000 

— - . 



... 


1 . 




„ . .. 
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TABLE II. - Concluded. VARIATION OF ROOT MEAN SQUARE MAGNITUDES OF MAGNETIC -FLUX AND 
ELECTRIC -FIELD STRENGTH WITH RADIAL POSITION 


[Radial position r normalized to 1 at boundary r = a.] 


J 

. Nondimensional 
radius, 
r/a 

i 

' Root mean 



Nondimensional radius to field penetration ratio, 

h 


square 

magnitudes 

0. 1 

0.5 

1.0 

1 2.0 

2. 5 

5.0 

7. 5 

10 

L ib _. 

20 

25 






a = » 1 







0 

/ 

0.998 

0.940 

0.789 

0.439 

[ 

0. 304 

0.037 

0.004 

-4 

4x10 

3xl0" 6 

2xlO~ 8 

2X10" 10 


A 

0 

0 

0 

0 

'° 

0 

0 

0 

0 

0 

0 

0. 1 

/ 

0.998 

0. 941 

0.791 

0. 443 

0. 309 

0.039 

0.004 

5xl0~ 4 

5xl0” 6 

5xl0~ 8 

6x10' 10 

; 

A 

. 100 

.097 

.089 

.063 

.050 

.011 

.002 

2xl0~ 4 

3x10® 

4xl0~ 8 

4x10" 10 

0. 2 

/ 

0.998 

0.943 

0.797 

0.456 

0.323 

0. 046 

0. 006 

8xl0' 4 

1Xl0 ’f 

3xl0' 7 

5x10' 5 



.200 

. 194 

. 178 

. 128 

. 102 

.023 

.004 

6xl0 -4 

lxlO -5 

2xl0~ 7 

4xl0 -9 

0.3 

/ 

0.998 

0.946 

0.807 

0. 479 

0. 348 

0.060 

0.010 

0.002 

5x10" 5 

2xl0' 6 

5x10' 6 


A 

. 300 

. 292 

. 268 

. 197 

. 160 

.040 

.008 

.001 

5x10“ 5 

lxlO' 6 

4x10' 8 

! 0.4 

/ 

0.998 

0.950 

0.821 

0. 512 

0. 384 

0. 083 

0. 018 

0. 004 

2xl0' 4 

lxlO" 5 

5xl0' 7 

J 

1 

A 

.400 

.390 

. 361 

. 272 

. 225 

.065 

.016 

.004 

2xl0' 4 

9x10' 6 

5x10' 7 

I 0. 5 

/ 

0.998 

0.955 

0. 839 

0. 555 

0.435 

0. 121 

0. 034 

0.010 

8xl0” 4 

-J 

X 

o 

1 

ai 

5x10' 6 


A 

1 .500 

.488 

.456 

. 355 

.300 

. 103 

.031 

.009 

8xl0 -4 

6xl0' 5 

5x10' 6 

I 

0.6 

/ : 

0.999 

0.962 

0.862 

0. 611 

0. 501 

0. 179 

0.065 

0.024 

0.003 

4xl0' 4 

6x10" 5 


A 

f .600 

. 588 

. 555 

.449 ! .390 

. . / ... 

. 162 

.062 ! 

.023 

.003 

4xl0' 4 

6X10" 5 

0. 7 

/ ! 

0.999 

0.969 

0.888 

0. 681 

0. 585 

0. 271 

0. 127 

0.060 

0.013 

0.003 

7xl0' 4 


A : 

.700 

. 689 

.658 

. 557 

.499 

. 255 

. 123 

.058 

.013 

.003 

7xl0' 4 

0. 8 

/ ; 

0.999 

co 

c- 

02 

d 

0.920 

0.768 

0. 693 

0.415 

0. 251 

0. 152 

0.056 

0.021 

0. 008 


/ 1 I .B00' 

.791 

.766 

. 682 

. 632 

.401 

. 246 

. 150 

.055 

L_ 

.020 

. 008 

0.9 

J I 

0. 999 1 

0.989 

0.957 

0.873 

0. 829 

0. 642 

0. 499 : 

0. 388 

0.235 

0. 143 

0. 086 


A 

. "j 

.900 

.895 

.880 

.828 { .796 

I_. , 

j . 632 

.495 

.386 

; 

. 235 

. 142 

.086 

1.0 

/ 

1. 000 

1.000 

1. 000 

1. 000 j 

1. 000 

1. 000 

1. 000 

1.000 

1.000 

1. 000 

1. 000 


A 

1. 000 , 

1 

t 

1.000 

1. 000 

1.000 1.000 

1. 000 

i 

1.000 

1.000 

1 

1.000 

i 

1.000 

1.000 
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TABLE m. - VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA AND THOSE 


AT BOUNDARY (r = a) a 


[ Nondimensional 

Phase 


Nondimensional radius to field penetration ratio, h 


radius, 

angle 










- 


r/a 


0. 1 

0.5 

1.0 

2.0 

2.5 

5.0 

7. 5 

10 

15 ; 20 

25 






a 

= 0 







0 

<P 

0. 145 

3.6 

14.3 

52.3 

74.7 

178.9 

280. 6 

382. 1 

564.9 

787. 5 

990. 2 


<p\ 

.072 

1.8 

7.2 

28.3 

43.4 

138. 6 

238. 6 

331.3 

556. 5 

743.6 

946.0 

0.1 

<P 

0. 142 

3. 5 

14. 1 

51.7 

73.8 

175.4 

272. 6 

367.9 

526. 3 

735.2 

915. 5 


<pi 

.071 

1.8 

7. 1 

28 

42.9 

136. 8 

234. 6 

332. 1 

525. 2 

715.3 

902. 6 

0. 2 

<P 

0. 138 

3.4 

13.7 

50 

71. 1 

164. 7 

249.5 

329.8 

468. 3 

649.3 

811. 2 



.069 

1.7 

6.9 

27. 1 

41. 6 

131.4 

222. 6 

311.0 

480.3 

644.7 

807.4 

0. 3 

<P 

0. 130 

3.3 

12.9 

47. 1 

66.6 

147.7 

210.0 

285.6 

426. 3 

567.9 

709. 5 


<p i 

.065 

1.6 

6. 5 

25.7 

39.4 

122.5 

203. 1 

278. 6 

422.6 

565. 1 

707.4 

0.4 

9 

0. 120 

3.0 

11.9 

43. 2 

60.4 

126.6 

184.0 

243.9 

365. 2 

486.6 

608. 1 


<p\ 

.060 

1.5 

6.0 

23.7 

36. 2 

110.3 

177.9 

240.4 

1362.8 

484.9 

606.7 

0.5 


0. 107 

2.7 

10.7 

38. 1 

52. 6 

104. 3 

152. 3 

203. 1 

304. 2 

405.4 

506.6 


<p\ 

.054 

1.3 

5.4 

21. 1 

32. 2 

95. 2 

149.4 

200.7 

302.7 

404.3 

505.8 

0.6 


0. 092 

2.3 

9. 1 

32.0 

43. 5 

82.4 

122.0 

162. 5 

223. 3 

324.3 

405. 3 


n 

.046 

1. 2 

4.4 

18.0 

27. 3 

77.8 

119.9 

160.8 

242. 3 

323. 6 

404.7 

0.7 

<p 

0. 073 

1.4 

7. 2 

24.9 

33. 2 

61. 3 

91.5 

121.8 

182. 5 

243.3 

303.9 



.037 

.9 

3.5 

14.3 

21. 6 

59. 1 

90. 1 

120.8 

181. 8 

242.7 

303. 6 

0. 8 

<p 

0.052 

1. 3 

5.1 

17. 1 

22.4 

40.7 

60.9 

81.2 

121. 6 

162. 1 

202. 6 


<pi 

.026 

.7 

2.4 

10.0 

15. 1 

39. 6 

60. 1 

80. 6 

121. 2 

161.9 

202. 5 

0.9 

V 

0.027 

0.7 

2.7 

8.8 

11.2 

20.3 

30. 5 

40.6 

60,8 

81. 1 

101. 3 


Vi 

.014 

.4 

1.4 

5.3 

7.9 

19. 9 

30. 1 

40. 3 

60. 6 

80.9 

101.2 

1.0 

<p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Vi 

o 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 


a FieIds in the plasma always lag behind those at r = a, and the absolute values of the phase 
angles in degrees are presented herein. 



TABLE ID. - Continued. VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA AND 


THOSE AT BOUNDARY (r = a ) a 


Nondimensional 

radius, 

r/a 

Phase 

angle 


Nondimensional radius to field penetration ratio, 

h 

! 

0. 1 

0. 5 

1.0 

2.0 

2. 5 

5.0 

7. 5 

1 10 

\ 

15 

j 20 

25 






a = 

1 







0 

9 

0. 102 

2. 50 

9. 30 

30. 5 

41. 6 

97.7 

152.8 

207.7 

317. 5 

427. 1 

530.0 



.051 

1. 26 

4. 99 

18.3 

27.0 

78.8 

135.2 

187. 5 

290.0 

419.0 

520.0 

0. 1 

9 

0. 100 

2.45 

9. 20 

29.7 

40.9 

95.2 

147.3 

198.4 

298.4 

397. 1 

495.2 



.050 

1. 24 

4. 89 

17.9 

26. 5 

76.7 

129.2 

181. 5 

275.0 

387. 1 

487.9 

0.2 

9 

0. 097 

2. 38 

8.90 

28. 5 

39. 1 

88.4 

133.8 

177.6 

254. 5 

351.6 

439. 1 


*i 

.049 

1.21 

4.72 

17.3 

25. 6 

73.0 

121. 3 

168.3 

259.4 

348. 5 

436.8 

0.3 

9 

0.092 

2. 25 

8.40 

26.6 

36. 1 

78. 7 

117.0 

154.7 

230.9 

307.4 

384.0 



.046 

1. 14 

4. 46 

16. 3 

24.0 

67. 2 

109.6 

150. 1 

228. 5 

307.8 

382.8 

0.4 

9 

0.085 

2. 07 

7. 70 

24.0 

32.2 

67. 6 

99.8 

132. 2 

197.7 

263.4 

329. 1 


n 

.043 

1.05 

4. 11 

15.0 

21.9 

59.8 

95.6 

129.7 

196.3 

262.4 

328.3 

0.5 

V 

0.076 

1.85 

| 6. 83 

20.8 

27.7 

56. 1 

82.9 

110.0 

164.7 

219.5 

274.2 


*1 

.038 

.94 

3. 66 

13.2 

19.3 

51. 1 

80.4 

i 

108. 5 

163.8 

218.8 

273.7 

0.6 


0.065 

1. 57 

5. 77 

17. 2 

22. 5 

44. 7 

66.2 

88.0 

131.7 

175. 5 

! 219.3 


*1 

.032 

.80 

3. 12 

11. 1 

16. 1 

41. 6 

64.7 

87.0 

131. 1 

175. 1 

219. 1 

0.7 

<p 

0.052 

1. 25 

4. 55 

13. 2 

17. 1 

33.4 

49.6 

65.9 

98.8 

131. 6 

164. 5 


*1 

.026 

.64 

2. 47 

8.8 

12. 6 

31.6 

48.7 

65.3 

98.4 

131.4 

164.3 

0.8 

<p 

0.036 

0. 88 

3. 17 

8.9 

11. 5 

22. 2 

33.0 

43.9 

65.8 

87.7 

109.7 


<p 1 

.018 

.45 

1. 74 

6. 1 

8. 7 

21. 2 

32.5 

43. 6 

65. 6 

87. 6 

109.5 

0.9 


0.019 

0.46 

1. 65 

4. 5 

5.7 

11. 1 

16. 5 

22.0 

32.9 

43.9 

54.8 


*1 

.010 

.24 

.91 

3.2 

4. 5 

10.7 

16.3 

21.8 

32.8 

43.8 

54.8 

1.0 

<p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


a Fields in the plasma always lag behind those at r = a, and the absolute values of the phase 
angles in degrees are presented herein. 



TABLE IB. - Continued. VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA AND 


THOSE AT BOUNDARY (r = a) a 


Nondi mensional 

Phase 


Nondi mensional radius to field penetration ratio 

h 


radius, 

angle 












r/a 


0. 1 

0. 5 

1.0 

2.0 

2,5 

5.0 

7.5 

10 

15 

20 

25 






Of = 2 







0 

<P 

0.064 

1.560 

5.81 

18.30 

25. 1 

58.7 

91.8 

124.8 

190.7 

256. 5 

320.0 



.064 

.790 

3. 10 

11. 40 

16.9 

47.8 

82.0 

112.0 

182.0 

260.0 

319.0 

0. 1 

<P 

0.063 

1. 540 

5.71 

18.00 

24.7 

57. 2 

88.4 

119.0 

179.0 

238.3 

297. 3 


9>1 

.032 

.785 

3.06 

11. 10 

16.3 

46.4 

77.8 

109. 2 

171. 2 

232.4 

292.9 

0.2 

<P 

0.061 

1.490 

5. 52 

17. 30 

23. 5 

53.0 

80. 2 

106. 5 

158.8 

211. 1 

263.6 



.031 

.761 

2.96 

10. 70 

15.7 

44. 1 

72.9 

101.0 

155.7 

209.2 

262. 3 

0. 3 

<P 

0.058 

1.410 

5. 21 

16. 00 

21.7 

47. 1 

70. 2 

92.9 

138.6 

184. 6 

230. 6 


<Pl 

.029 

.721 

2.80 

10. 10 

14.7 

40. 5 

65.8 

90.0 

137. 1 

183. 6 

229.9 

0.4 

<P 

0.054 

1. 300 

4.77 

14. 40 

19.3 

40. 5 

59.9 

79.4 

118.7 

158. 1 

197. 6 



.027 

. 665 

2. 58 

9. 23 

13.4 

35.9 

57.3 

77.8 

117.8 

157. 5 

197. 1 

0. 5 


0.048 

1. 160 

4. 22 

12. 50 

16. 5 

33.7 

49.8 

66. 1 

98.9 

131.7 

164.6 



.024 

.594 

2. 30 

8. 14 

11.8 

30. 7 

48. 2 

65. 1 

98. 3 

131.3 

164.3 

0. 6 


0.041 

0. 987 

3. 56 

10. 30 

13.4 

26.8 

39. 7 

52. 8 

79. 1 

105.3 

131.7 


Vi 

.020 

. 506 

1.95 

6.85 

9.8 

24.9 

38. 8 

52. 2 

78.7 

105. 1 

131. 5 

0.7 

<p 

0.033 

0.784 

2. 80 

7.87 

10. 2 

20.0 

29. 8 

39. 6 

59.3 

79.0 

98.8 


^1 

.016 

.403 

1. 55 

5. 37 

7. 6 

18.9 

29. 2 

39. 2 

59. 1 

78.9 

98.6 

0.8 

(P 

0.023 

0. 551 

1.95 

5.33 

6.8 

13.3 

19.8 

26.4 

39. 5 

52.7 

65.8 


<P1 

.012 

. 284 

1.09 

3.72 

5. 3 

12.7 

19. 5 

26. 2 

39.4 

52. 6 

65.8 

0.9 


0.012 

0. 290 

1.01 

2. 69 

3.4 

6.7 

9.9 

13. 2 

19. 8 

26.3 

32.9 


<pl 

.006 

. 150 

. 57 

1.92 

2.7 



6.4 

9.8 

13. 1 

19. 7 

26. 3 

32.9 

1.0 

<p 

0 

0 

0 

0 

0 

0 

i 0 

0 

0 

0 

0 



0 

0 

0 

0 

. 

1° 

0 

0 

0 

i 

0 

0 

0 


a Fields in the plasma always lag behind those at r = a, and the absolute values of the phase 
angles in degrees are presented herein. 
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TABLE III. - Continued. VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA 


AND THOSE AT BOUNDARY (r = a) a 


Nondime ns ional 
radius, 
r/a 

Phase 

angle 


Nondimensional radius to field penetration ratio, h 


0. 1 

0.5 

! 1.0 

2.0 

i 

j 2. 5 

5.0 

7. 5 

10 

15 

20 

25 






a = 3 








0 

<P 

0. 047 

1. 120 

4. 08 

12. 90 

17. 5 

41.0 

64.0 

87.0 

133.0 

178.9 

228. 1 


V 1 

.022 

.564 

2. 20 

7.90 

11. 6 

33. 2 

56.7 

79. 6 

121. 5 

170.0 

221.0 

0. 1 

<p 

0.045 

1.090 

4.02 

12. 60 

17.2 

39.9 

61.7 

83.0 

124.8 

166. 1 

207. 3 


<p 1 

.022 

.555 

2. 16 

7. 80 

11.4 

32.4 

54.3 

76. 2 

119.4 

162.0 

206. 2 

0.2 


0. 043 

1.050 

3. 88 

12. 10 

16.4 

36.9 

55.9 

74.3 

110. 7 

147.2 

183.8 



.022 

. 538 

2. 09 

7. 53 

11.0 

30.8 

50.8 

70.4 

108. 5 

145.8 

182.8 

0. 3 

<p 

0.041 

0.998 

3.66 

11. 20 

15. 1 

32.8 

48.9 

64.8 

96. 6 

128.7 

160.7 


V 1 

.021 

. 510 

1. 98 

7.09 

10.3 

28.3 

45. 9 

62.8 

95. 6 

128.0 

160. 2 

0.4 

<p 

0.038 

0.920 

3.35 

10. 10 

13. 5 

28. 2 

41. 8 

55. 3 

82.8 

110. 2 

137. 7 


<pl 

.019 

.470 

1.82 

6.48 

9.4 

25. 1 

40.0 

54. 2 

82. 1 

109.8 

137.4 

0.5 

cp 

0.034 

0.820 

2.96 

8.71 

11. 5 

23.4 

34.7 

46. 1 

68.9 

91. 8 

114.8 


n 

.017 

.420 

1. 62 

5.71 

8. 2 

21.4 

33. 6 

45.4 

68. 5 

91.6 

114.6 

0.6 

v 

0.029 

0. 700 

2. 50 

7. 16 

9.4 

18.7 

27.7 

36.8 

55. 1 

73.5 

91.8 


<?! 

.014 

.358 

1. 38 

4.80 

6,9 

17.4 

27.0 

36.4 

54.9 

73. 3 

91.7 

0.7 

<P 

0.023 

0. 553 

1.96 

5.48 

7, 1 

14.0 

20.8 

27. 6 

i 41.3 

55, 1 

68.9 


*1 

.012 

. 285 

1.09 

3.76 

5.3 

13. 2 

20.4 

27. 3 

41.2 

55.0 

68.8 

0.8 

V 

0.016 

0. 389 

1. 37 

3.71 

4.8 

9.3 

13.8 

18.4 

27. 6 

36.7 

45.9 


<p 1 

.008 

. 201 

.77 

2.61 

3.7 

8.9 

13. 6 

18.2 

27.5 

36.7 

45. 6 

0.9 

<P 

0. 009 

0. 204 

0.71 

1.87 

2.4 

4.6 

6.9 

9. 2 

13.8 

18.4 

22.9 


<P 1 

.004 

. 106 

.40 

1.35 

1.9 

4. 5 

6.8 

9. 1 

13.7 

18.3 

22.9 

1.0 

<p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


*1 

0 j 

0 

0 

0 

0 

0 

0 

0 

° l 

0 



a Fields in the plasma always lag behind those at r = a, and the absolute values of the phase 
angles in degrees are presented herein. 



TABLE HI. - Continued, 


VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA 


AND THOSE AT BOUNDARY (r = a) a 


Nondimensional 

radius, 

r/a 

Phase 


Nondimensional radius to field penetration ratio, h 


angle 

0. 1 

0. 5 

1.0 

2.0 

2.5 

5.0 

7. 5 

10 

15 

20 

25 






a = 5 








0 

<P 

0.029 

0. 685 

2. 52 

7. 87 

10. 80 

25.2 

39.4 

53.4 

81.8 

110. 1 

138.0 


<p i 

.014 

.352 

1. 37 

4.91 

7. 17 

20. 5 

35.4 

49.6 

78. 3 

109. 7 

136.0 

0. 1 

<P 

0.028 

0. 675 

2. 48 

7,76 

10.62 

24. 5 

37.9 

51.0 

76.8 

102. 2 

127. 5 



.014 

.344 

1. 34 

4.82 

7.06 

19.9 

33.4 

46.9 

73. 5 

99.7 

125.6 

0. 2 

<P 

0.027 

0. 654 

2. 40 

7. 43 

10. 10 

22. 7 

34.4 

45.7 

68. 1 

90. 6 

113. 1 


*i 

.014 

.334 

1. 29 

4.65 

6.80 

18. 9 

31.3 

43.. 3 

66. -8 

89. 7 

112. 5 

0. 3 

<P 

0.026 

0. 619 

2. 26 

6. 90 

9.39 

20. 2 

30. 1 

39.8 

59.5 

79. 1 

98.9 



.013 

. 316 

1. 22 

4. 38 

6. 37 

17.4 

28. 2 

38. 6 

58.8 

78. 7 

98. 6 

0.4 

<P 

0.024 

0. 570 

2. 07 

6. 20 

8. 28 

17. 3 

25.7 

34. 1 

50.9 

67. 8 

84.7 


<pi 

.012 

. 292 

1. 13 

4.00 

5.79 

15.4 

24. 6 

33.4 

50. 5 

67.6 

84.5 

0. 5 

V 

0.021 

0. 507 

1.83 

5. 36 

7.07 

14.4 

21. 3 

28.3 

42.4 

56. 5 

70. 6 


*1 

.011 

. 260 

1.00 

3. 53 

5.08 

13. 1 

20. 7 

27.9 

42. 2 

56.3 

70. 5 

0.6 

(P 

0.018 

0.432 

1. 54 

4.40 

5.75 

11. 5 

17.0 

22.7 

33.9 

45. 2 

56. 5 


Vi 

.009 

. 222 

.85 

y 2. 96 

4. 24 

10. 7 

16.6 

22.4 

33. 8 

45. 1 

56.4 

0.7 

<p 

0.014 

0. 343 

1. 21 

3. 37 

4.36 

8. 6 

12. 8 

17.0 

25.4 

33. 9 

42.4 


*1 

.007 

. 176 

.68 

2. 32 

3.30 

8. 1 

12. 5 

16.8 

25. 3 

33.8 

42.3 

0.8 

<p 

0.010 

0. 241 

0. 84 

2. 28 

2.92 

5.7 

8.5 

11. 3 

17.0 

22.6 

28. 2 


*1 

.005 

. 124 

.47 

1. 61 

2. 27 

5. 5 

8.4 

11. 2, 

16.9 

22.6 

28.2 

0.9 

<P 

0.005 

0. 126 

0. 44 

1. 15 

1.46 

2.9 

4.3 

5. 7t 

8. 5 

11. 3 

14.1 


V 1 

.003 

.066 

.25 

.83 

1. 16 

2.7 

4. 2 

5.6 

8. 5 

11. 3 

24. 1 

1.0 

<p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


a Fields in the plasma always lag behind those at r = a, and the absolute values of the phase 
angles in degrees are presented herein. 



TABLE HI. - Continued. VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA 


AND THOSE AT BOUNDARY (r = a) a 


Nondimensional 

Phase 


Nondimensional radius to field penetration ratio, 

h 


radius, 

r/a 

angle 

0. 1 

0. 5 

1.0 

2.0 

[ 2.5 

5.0 

I ... 7 * 5 

L io _ 

1 15 

20 

25 






a = 10 








0 

<P 

0.014 

0. 349 

1. 280 

4.02 

5.46 

12.70 

19.90 

27.10 

41.4 

55.6 

71.8 


*1 

.007 

.178 

.681 

2.49 

3.65 

10. 40 

17. 50 

24.50 

38.7 

54.3 

70.0 

0. 1 

<P 

0.014 

0. 342 

1. 260 

3.93 

5.37 

12. 40 

19. 20 

25.80 

38.8 

51.7 

64. 5 


<pi 

.007 

. 175 

.678 

2.44 

3. 57 

10. 10 

16. 90 

23.70 

37. 1 

50.4 

63.5 

0.2 

<P 

0.014 

0.331 

1. 220 

3.76 

5. 11 

11. 50 

17. 40 

23. 10 

34.4 

45.8 

57.2 


<p i 

.007 

. 169 

.656 

2.36 

3.44 

9.58 

15. 80 

21.90 

33.7 

45.3 

56.9 

0. 3 

<P 

0.013 

0.314 

1. 150 

3.49 

4.71 

10. 20 

15. 20 

20. 10 

30. 1 

40.0 

49.9 


<p i 

.007 

. 160 

.621 

2. 22 

3. 23 

8.79 

14. 30 

19. 50 

29.7 

39.8 

49.8 

0.4 

<P 

0.012 

0.289 

1.050 

3. 14 

4.18 

8.76 

13.00 

17. 20 

25.7 

34.3 

42.8 


*1 

.006 

. 148 

.572 

2.03 

2.93 

7. 79 

12. 40 

16.90 

25.5 

34.2 

42.7 

0.5 

9 

0.011 

0. 257 

0.927 

2.71 

3. 58 

7. 30 

10. 80 

14. 30 

21.4 

28.6 

35.7 


<Pl 

.005 

. 132 

. 508 

1.79 

2. 57 

6.65 

10. 50 

14. 10 

21.3 

28. 5 

35.6 

0.6 

9 

0.009 

0. 219 

0.782 

2. 23 

2.91 

5.81 

8. 62 

11. 50 

17. 1 

22.8 

28.6 


Vl 

.005 

. 112 

.432 

1.50 

2. 14 

5. 40 

8.41 

11. 30 

17. 1 

22.8 

28. 5 

0.7 

9 

0.007 

0. 174 

0.614 

1.70 

2. 20 

4.35 

6.45 

8.58 

12.9 

17. 1 

21.4 


[ *1 

.004 

.089 

.343 

1. 18 

1.67 

4. 10 

6.33 

8.50 

12.8 

17. 1 

21.4 

0.8 

9 

0.005 

0. 122 

0.427 

1. 15 

1.48 

2. 89 

4. 30 

5.72 

8. 6 

11.4 

14.3 


9 X 

.003 

.063 

. 240 

.81 

1.15 

2.75 

4. 23 

5. 67 

8. 5 

11.4 

14. 3 

0.9 

9 

0.003 

0.064 

0. 221 

0. 58 

0.74 

1.44 

2. 15 

2. 86 

4.3 

5.7 

7. 1 


V 1 

.001 

.033 

. 126 

.42 

.59 

1. 38 

2. 12 

2. 84 

4.3 

5.7 

7. 1 

1.0 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 


a Fields in the plasma always lag behind those at r = a, and the absolute values of the phase 
angles in degrees are presented herein. 



TABLE m. - Continued. VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA 


AND THOSE AT BOUNDARY (r = a) a 
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TABLE HI. - Concluded. VARIATION OF PHASE ANGLES BETWEEN FIELDS IN PLASMA 


AND THOSE AT BOUNDARY (r = a) d 


N ondi mens ional 

Phase 


Nondimensional radius to field penetration ratio, 

h 


radius. 

angle 

— 

— 

— 


r 

— 

— 


- ■■ -- 

T 


r/a 

0. 1 

0.5 

1.0 

a 

2.0 
= 100 

2. 5 

5.0 

7. 5 

j 10 

15 

j 20 

25 

0 

<P 

0.0010 

0. 035 

0. 128 

! 

0. 401 

0. 548 

! 1. 28 

r 

1.99 

2.72 

4. 15 

5. 58 

7. 19 


*i 

.0007 

.019 

.069 

! .250 

1 

. 370 

1.04 

1 

1.78 

2. 40 

4. 00 

5. 50 

7.06 

0. 1 

<P 

0.0010 

0.034 

0. 126 

i ; 

: 0. 394' 0. 539 

j 1.24 

1.92 

2. 59 

3.90 

5. 18 

6. 47 



.0007 

.018 

.068 

, .245 

• .360 

i 

] 1.01 

1.70 

2. 38 

3.73 

5.06 

6. 37 

0.2 

<P 

0. 0010 

0.033 

0. 122 

0. 377 

0. 513 

j 1. 15 

1. 74 

2.32 

3.46 

4. 59 

5.74 


V 1 

.0007 

.017 

.066 

. 237 

. 350 

.96 

1.59 

2. 20 

3. 39 

4. 55 

5.71 

0.3 

<p 

0. 0010 ; 0. 032 

0. 115 

0. 350 

0. 472 

! 1.02 

1. 51 

2.02 

3.02 

4.02 

5.02 


*1 

.0007 

.016 

.062 

. 223 

. 320 

.88 

i 

1.43 

1.96 

2.98 

3.99 

5.00 

0.4 

cp 0.0010 

0.029 

0. 105 i o. 314 

0. 420 

0.88 

1. 30 

1.73 

2. 58 

3. 44 

4.30 


<Px ; 

.0006 

.015 

.057 

. 203 

. 290 

.78 

1.25 

1. 69 

2.56 

3.43 

4. 29 

0. 5 

<P ! 

0. 0010 

0.026 

0.093 

0. 272 

0. 359 

0.73 

1.08 

1.44 

2. 15 

2. 87 

3. 58 


: 

.0005 

.013 

.051 

. 179 

. 260 

. 67 

1.05 

1.42 

2. 14 

2.86 

3. 58 

0. 6 

<P 

0.0010 

0.022 

0.079 

0. 220 

0. 292 

0. 58 

0. 86 

1. 15 

1.72 

2. 29 

2.87 


<?! , 

.0005 

.011 

.043 

. 151 

. 210 

. 54 

.84 

1. 14 

1.71 i 

2. 29 

2.86 

0.7 

<P 1 

0. 0010 

0.017 

0.062 

0. 170 

0. 221 

0.44 

0. 65 

0. 86 

1. 29 

1.72 

2. 15 


Vi ■ 

.0004 

.009 

.034 

. 118 

. 170 

.41 

. 64 

.85 

1. 29 

1.72 

2. 15 

0.8 

<p 

0. 0010 

0.012 

0. 043 

0. 120 

0. 148 

!0. 29 

0.43 

0. 57 

0.86 

1. 15 

1.43 


<P1 

.0003 

.006 

.024 

.082 

. 120 

.28 

.42 

. 57 

.86 

1. 14 

1.43 

0.9 

<p 

0. 0003 

0. 006 

0, 022 

0. 060 

0. 074 

0. 14 

0. 22 

0. 29 

0. 43 

0. 57 

0.72 


<Px 

.0004 

.003 

.013 

.042 

.060 

. 14 

.21 

. 28 

.43 

. 57 

.72 

1.0 

<p 

0 

0 1 

0 

0 

0 

0 

0 

0 

0 

0 ,0 


<P\ 

0 

i 

0 jo 

0 

0 I 

0 

0 

0 

0 

0 jo 


^Fields in the plasma always lag behind those at r = a, and the absolute values of the phase 
angles in degrees are presented herein. 










APPENDIX C 


SOLUTIONS WHEN a»l 
Solutions for Nondimensional Power 

Using equations (30) and (31), the nondimensional power given by equation (34) may 
be rewritten as 


g _ 6 (l + Q? 2 ) 


1/4 


Re 


4a 


sin 


'cot -1 ® 


/-icot-V 

exp ( — — 


J o({a> 


(Cl) 


When a » 1, the following relations hold 


i-l. 


expf- 1 cot ~ XQ? V^exp^ 


sin/^ cot -— 1 


V2 

_i_ 

2a 


! a = i exp/— ► hi expf—\ 

a 6 [2aJ K2a) 


(C2) 


(C3) 


(C4) 


Substituting these values into equation (Cl) yields 


P = 


a 


Re 


r 

J i 

hi expf—) 

\2a/J 

J 0 

hi exp/——) 
\2a/_ 


ex *(-k) 

j 


Using the series definition of the Bessel functions yields, when a » 1, 

J »[ hlexp &)] e ‘ I ° (h, + ^ I ' <h) 


and 


r i[ hi exp (£>] “ exp &) [ 1Ii<h) ■ ^ l2<h) ] 


(C5) 


(C6) 


(C7) 
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Equation (C5) becomes, therefore 


P = 


a 


Re. 


yft h 


i 0 (h) + Jii ij(h) 

u 2a 1 


exp 


te) 


(C8) 


where the I(h)'s are the modified Bessel functions. Reducing equation (C8) and rational- 
izing yield 


P = 


2V2 


q(h) - I Q (h)I 2 (h) 


Iq(h) + ~ I 2 (h) 
4 a 2 


(C9) 


For a » 1, the second term in the denominator will be very small compared with 
2 

IqOi) and equation (C9) reduces to 


P = 


2V2 


Ij(h) I 2 (h) 


(38) 


Solution for Change in Inductance 


The solution for AL/L is given by equations (45) and (42) as 


AL 


op 


= 1 - ®/f(l ♦ al 1/2 lm[^l - [(l + a , 2 ) 

a \ L J L 


1/2 - 11/2 


arj Re 


W 


Wr 


Using equations (35) and (18) and simplifying yield 


(Cll) 


Ak. 


op 


W2 cos/ cot 'Mlm 

piM 

- 'ft sh/ cot " lQ! W 

" J l<y]l 

1 V 2 / 

_V«a>J 

w ' \ 2 / 

_ J 0<«a>Jj 


(C12) 
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or 


AL * 2 * 

=1 — Im 


op 


L J 0<«a> 


exp^- - cot 


For a » 1, using equations (C2) and (C4) yields 

exp fe)] 


AL 4 2 T 

— — = 1 - - Im 

L h 

op 


r o[ hi exp (^). 


exp 


(-=) 


Using equations (C6) and (C7) and rationalizing yields 


IoWl^h) + — I^hJIgOi) 

AL _ j _ 2 4 oft 


op 


l 2 0 (h) + I^h) 

4ctT 


or for a » 1 


AL = L_ 


op 


! _ 2 r i( h ) 


L h V h >J 


(C 13) 


(C14) 


(C 15) 


(46) 
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APPENDIX D 


DETERMINATION OF CIRCUIT PARAMETERS AND CIRCUIT CHECK 


The measurement of the coil current and inductance was complicated by the fact that 
the operating frequency of the system is near the self -resonant frequency of the coil when 
the stray capacitance of the Faraday shield is taken into consideration, hi order to de- 
termine whether or not a true measurement of the coil current and inductance was ob- 
tained, it was necessary to determine the values of the individual circuit parameters. 
With this knowledge a check may be made of the validity of the various current and volt- 
age measurements by checking the system when no plasma is present. The circuit pa- 
rameters that must be determined are the coil inductance Lq, the lead inductance from 
the coil to ground, the inductance of the strap connecting the variable capacitors plus the 
lead to the coil inductance Lg, and the stray capacitance C g (fig. 5, p. 16). 

The inductance of the coil was calculated by using the equation (Welsby, ref. 13) 


47r 2 N 2 b 2 -9 

K n C c xlO y microhenries 


J 0 


l 


(Dl) 


where K is a correction factor used for short coils 
n 





(D2) 


and C c is a correction factor used for the spacing of the coil turns 


C 


c 


1 _ £(A’ + B’) 
?rbNK n 


(D3) 


where A' = 2. 3 log^d'/c, B' = 0. 366[l - 2. 5/N + 3. 8/N 2 ], d’ is the diameter of the 
coil tubing, and c is the winding pitch of the coil. 

The calculated inductance of the coil using the mean coil radius is Lq = 1.0 micro- 
henries. The coil inductance was measured at 17. 5 megahertz using a Q meter and an 
RX meter. The measured Lq values were 1. 08 and 1. 05 microhenries, respectively. 

The inductances of the remaining portions of the circuit were obtained from mea- 
surements made using an L-C meter and a low -inductance bridge. The stray capacitance 
of the coil and Faraday shield was obtained by adding capacitors of known value across 
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I 


l 2 A rms 



Figure 12. - Circuit parameters and measurements. 

2 

the coil and measuring the resonant frequency c o . If 1/w is plotted against the added 
capacitance, the stray capacitance may be obtained from the intercept on the abscissa 
axis and the inductance is the slope of the curve. The measured values of the circuit pa- 
rameters as well as an indication of the measurements made during the experiments and 
tare data check are shown in figure 12. 

When the system is tuned for zero reflected power at 17. 5 megahertz the following 
relations hold for this circuit 


C s = 


6. 48 V (kV) 

L = ±- microhenries 

V A rms^ 

12. 74 [lg(A rms ) " V A rms^] 


V p (kV) 


picofarads 


J eff 


VfkV) 3 

— E xlO J ohms 


^ V A rms^ 


(D4) 

(D5) 

(D6) 


Where Z e ^> is the effective impedance of the L-R-C g portion of the circuit and was 
given by equation (48) as 


I 
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J eff 


coL 


ap 


(D7) 


If one defines the effective inductance of the circuit as = L a ^ + Lg then 

is (from eq. (52)) 


2. 5X1Q - 3 Cl (pF) | 82 . 78 
3.02X10 -5 Cj(pF) + 1 C 2^ pF ^ 


microhenries 


(D8) 


The system is checked by taking a series of measurements with no plasma present. 
These tare results are presented in table IV. It can be seen from this table that the 


TABLE IV. - CIRCUIT CHECK DATA WHEN NO PLASMA IS PRESENT 


Root mean 

Root mean 

Peak 

Inductance, 

Stray 

Calibrated vari- 

Effective 

Effective 

Apparent 

' Inductance 

square 

square coil 

voltage, 

L, 

capacitance, 

able capacitator 

inductance, 

impedance, 

inductance, 

1 of strap 

current, 

% 

A 

current, 

I 3- 

A 

v p , 

V 

m h 

C s’ 

PF 

c i> 

pF 

c 2 , 

pF 

L eff> 

mh 

Z efF 

ohms 

L aP’ 

pH 

plus lead 
inductance 
of coil, 

rms 

rms 







1 

I 

! 




(a) 

(b) 



(c) 

(d) j 

(e) 

pH 

7 

12.5 

2165 

1. 122 

32. 33 

762 j 

32.6 

I 2.64 

221,8 

2.02 

j 0.62 


8 

9 

10 

11 

12 

13 

14 

15 

Average 


14.5 
16. 3 
18.3 
20 
22 

23.8 
25. 5 
27. 5 


2510 

2810 

3190 

3495 

3815 

4110 

4405 

4760 


Measured value 


1.121 
1.117 
1.129 
1. 132 
1.123 
1.119 
1. 19 
1 . 121 

1. 122 

- 1 . 1 


32.96 
33.07 
33. 10 
32. 70 
33.30 
33.40 
33.20 
33.40 

33.05 

~32 


224 
223 
228 
227 
227 
226 

225 
227 

225 


2. 04 
2. 03 
2. 07 
2. 07 
2. 07 
2. 06 
2.05 
2. 06 

2.05 


.60 
.61 
.57 
.57 
. 57 
. 58 
.59 
.58 

0.588 

0.58 


a See eq. (D4). 
b See eq. (D5). 
c See eq. (D8). 
d See eq. (D6). 
e See eq. (D7). 


values of the circuit parameters obtained from the tare data agree quite well with the 
measured values. Therefore, we believe that the measured coil current and inductance 
are accurate. 

In order to obtain the data when the plasma is present, the current Ig is set and 
the power coupled to the plasma is obtained. The measurement of the change of induc- 
tance with a plasma is difficult because of the small magnitude of this change. It was 
for this reason that the change in coil inductance was determined using two methods. The 
first method (eq. (47) in text) used the V and Ig measurements. The one shortcoming 
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of this method is that the L value calculated from the tare data varied slightly from day 
to day for various coil currents. In an attempt to overcome this problem, series of data 
were taken at constant coil current. The measured values of the series and shunt capaci- 
tors did remain constant for the daily tare data check when the system was tuned for zero 
reflected power. The second method used to calculate the inductance change was to cal- 
culate from these capacitor measurements (see eq. (52) in text) and obtain AL 

from equation (50). 
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APPENDIX E 


DETERMINATION OF PLASMA RESISTIVITY 


The plasma resistivity is obtained by assuming that the resistivities due to electron- 
atom and electron-ion collisions are additive; that is, 

^^ei + ^ea < E1 ) 

When all ions are assumed to be singly charged, Spitzer (ref. 6) gives an equation for 
rj e ^ which becomes in SI units 

T] . = 65. 3 — ohm -meter (E2) 

61 t 3/2 

e 

where A = 127r(eQkT e )^ // ^/q^n^^ and is tabulated in reference 6. 

The resistivity due to electron-atom collisions is obtained from the equation 


ri 


ea 


mu' „ 
e ea 



(E3) 


where v' is the normalized electron -atom momentum -transfer collision frequency and 
is defined as (see ref. 14 or 15) 



where f . is the fraction ionized defined as n g /nQ. 

In equation (E4), the momentum -transfer cross section for electron-atom collisions 

was taken from Brown (ref. 16) as a function of the kinetic energy of the electrons. The 

resistivity r\ is plotted in figure 13 as a function of electron kinetic temperature with 

the fraction ionized as a parameter. In this plot, the In A was calculated using 
17 

n g = 1x10 per cubic meter in the rj e ^ calculation, an assumption which causes little 
effect since In A varies slowly with n g . 
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Figure 13. - Plasma resistivity as function of electron temperature 
with fraction ionized as parameter. 


In order to determine the resistivity for the plasma under consideration, the elec- 
tron density, the electron kinetic temperature, and the neutral density must be known. 

The electron density was obtained from the microwave interferometer data. 

The electron kinetic temperature was obtained from the relative intensities of helium 
singlet and triplet spectral lines using the technique presented in reference 12 for tenuous 
helium plasmas. The 5047 and 4713 A (504. 7 and 471. 3 nm) spectral lines were used, 

and the excitation coefficients for these lines were taken from reference 17. These spec- 

17 

tral lines will give accurate results for electron densities below 3x10 per cubic meter 
and will exhibit a minimum error for 3x10 17 < n g < 10 19 per cubic meter (see Podgornyi, 
ref. 18). A correction was made for cascading effects, and a study of the light emitted 
from a number of spectral lines indicated that the plasma did indeed fall into the tenuous 
regime (ref. 19). The relative wave length response of the entire system was calibrated 
using a standard ribbon filament tungsten lamp. 

The neutral density was obtained from the McLeod gage measurements. Since the 
pressures in the McLeod gage and the apparatus are not necessarily equal tinder the near 
free molecular flow conditions which were present, the expression used by Chubb and 
Seikel (ref. 14) 
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was used to calculate riQ from the measured pressure. In equation (E5), Pq is the 
McLeod gage pressure (N/m^), Tq and Tq are neutral temperatures in the McLeod 
gage and discharge. Neutral temperature Tq, which should be nearly room tempera- 
ture, was chosen to be 300° K. The neutral temperature of the discharge was assumed 
to be 450° K. 
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